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Abstract 

We analyze, within the wavelet theory framework, the wandering over a screen 
of the centroid of a laser beam after it has propagated through a time-changing 
laboratory-generated turbulence. Following a previous work (Fractals 12 (2004) 223) 
two quantifiers are used, the Hurst parameter, H, and the Normalized Total Wavelet 
Entropy, NTWS. The temporal evolution of both quantifiers, obtained from the laser 
spot data stream is studied and compared. This allows us to extract information of 
the stochastic process associated to the turbulence dynamics. 
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1 INTRODUCTION 



The purpose of this work is to statistically describe laser beam propagation 
through time-changing laboratory- generated turbulence. To do so, we ana- 
lyze data stream corresponding to the centroid position of the laser spot by 
using two different quantifiers obtained from the wavelet theory: the Hurst 
parameter, H, and the Normalized Total Wavelet Entropy, NTWS. The for- 
mer quantifier results from modeling the centroid's coordinates as a fractional 
Brownian motion (fBm) at stationary turbulence strength [1], while the latter 
has been used for a wider set of stochastic processes — see Ref. [2]. 

The fBm was discovered by Kolmogorov [3] and defined by Mandelbrot and 
Van Ness [4] as the only one family of processes which are gaussian, self-similar, 
and with stationary increments. The normalized family of these gaussian pro- 
cesses, , is the one with B^{0) = almost surely, E[i?'^(t)] = 0, and 
covariance 



E 



s\ 



\2H 



for s, t G M. The power exponent H is also known as scaling exponent and its 
range is bounded between and 1. The estimation of this parameter plays a 
key role modeling a fBm time series. One remarkable property of this family 
B^ is that the H parameter regulates the presence or absence of memory. 
In fact, it can be separated in three subfamilies accordingly: long-memory for 
1/2 < H < 1, memoryless at H = 1/2 (ordinary Brownian motion), and 
short-memory in the case < H < 1/2. Likewise, the Hurst parameter can 
be thought as the probability of the next increment of the signal have the 
same sign as the previous increment. Thus, it tunes the trajectory regularity. 
Fractional Brownian motions are continuous but non-differentiable processes 
(in the usual sense), and only give generalized spectra 1//" with exponents 
a between 1 and 3. As a nonstationary process, the fBm does not have a 
spectrum defined in the usual sense; however, it is possible to define a power 
spectrum of the form [5] : 

SbhU) = JJ^H+T- (2) 

Remember that this equation is not a valid power spectrum in the theory of 
stationary processes since it is a non-integrable function in the classical sense. 

Several properties evidence that wavelet analysis is well-suited to fBm: 

(1) fBm is nonstationary but the wavelet coefficients form a stationary pro- 
cess at each scale [6, 7]; 

(2) fBm exhibits positive long-range correlation in the range 1/2 < < 1 
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but wavelet coefficients have a correlation which is highly small as soon 
as > if + 1/2, where is the number of vanishing moments associated 
to the mother wavelet ip(t) [7, 8]; 
(3) the self-similarity of fBm is reproduced in its wavelet coefficients, whose 
variance varies as a power law as a function of scale j [6, 7] 



These features are due to the fact that the wavelet family tpafi is generated 
by dilations and translations of a unique admissible mother wavelet iplt). So, 
the family itself exhibits scale invariance. It should be noted that the first two 
properties are valid for any process with stationary increments [9]. 

In particular these properties are widely used for estimating H or the related 
spectral exponent a = 2H + 1 [5, 10, 11]. Through the Logscale Diagram the 
threefold objective: detection, identification and measurements of the scaling 
exponent can be achieved [12]. Basically the estimation problem turns into a 
linear regression slope estimation. 

In order to model the time-changing turbulence situation we consider a gener- 
alization where the parameter H is no longer constant, but a continuous func- 
tion of the time t [H — > H{t)). This generalization was introduced in financial 
research to model the behaviour of stock market index time series [13-15]. 
Also, it was recently used to characterize dynamic speckle or biospeckle [16]. 
A single scaling exponent would be unable to show the complex dynamics 
inherent to the data. The constraint of stationary increments is relaxed in this 
case. Multifractional Brownian motion (mBm) [17] was formalized as a class 
of processes which satisfies these properties. 

We calculate the time- dependent Hurst exponent by using the wavelet prop- 
erties. Provided that variations of H are smooth enough, the signal is divided 
into i non-overlapping temporal windows and the scaling exponent is cal- 
culated for each subset according to the procedure described at Ref. [1]. A 
sequence of Hurst parameter values is obtained. They give the local scaling 
exponent around a given instant of time. Artificially mBm were analyzed in 
order to test the quality of our estimator. In Figure 1 one can compare the 
theoretical and experimental results for a mBm with H changing linearly from 
0.1 to 0.9 with t. The Matlab code introduced by J. F. Coeurjolly [18] was 
implemented to simulate the mBm. The signal was divided in 64 temporal 
windows of 512 data points. We used the orthogonal cubic spline functions as 
mother wavelet and the resolution levels from j = — 9 to j = — 1. 

At the same time, the NTWS is also applied to study this time-changing 
turbulence. Introduced as a measure of the degree of order-disorder of the 
signal [2], it provides information about the underlying dynamical process 




(3) 
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associated with the signal. We define the NTWS as 

SwT = - E (4) 

j=-N 

where 

S'^^- = \nN (5) 
with A^, at least, the base 2 logarithm of the number of data points and {pj} 
represent the Relative Wavelet Energy (RWE). These are defined as 

Pj = Sj/Stot, (6) 

with £j = K[Cj{k)] the energy at each resolution levels j = —N, ■ ■ ■ , —2, — 1 
and Stot = J2j<o^j- They yield, at different scales, the probability distribution 
for the energy. It should be remarked that an orthogonal mother wavelet must 
be used within this theory — further details can be found at Ref. [1]. Indeed, 
a very ordered process can be represented by a signal with a narrow band 
spectrum. A wavelet representation of such a signal will be resolved in few 
wavelet resolution levels, i. e., all RWE will be (almost) zero except at the 
wavelet resolution levels which includes the representative signal frequency. 
For this special levels the RWE will be almost equal to one. As a consequence, 
the NTWS will acquire a very small, vanishing value. A signal generated by a 
totally random process or chaotic one can be taken as representative of a very 
disordered behavior. This kind of signal will have a wavelet representation 
with significant contributions coming from all frequency bands. Moreover, one 
could expect that all contributions will be of the same order. Consequently, 
the RWE will be almost equal at all resolutions levels, and the NTWS will 
acquire its maximum possible value. Higher values for wavelet entropy means 
higher dynamical complexity, higher irregular behaviour and, of course, lower 
predictability. 

The time evolution of NTWS can be easily implemented. So, it is widely used 
to study a wide set of nonstationary natural signal. In particular, it was in- 
troduced to quantify the degree of disorder in the electroencephalographic 
epileptic records giving information about the underlying dynamical process 
in the brain [19], more specifically of the synchrony of the group cells involved 
in the different neural responses. Also, monthly time series of different solar 
activity indices (sunspot numbers, sunspot areas and flare index) were ana- 
lyzed [20, 21]. The disorder content of solar cycle activity can be derived by 
analyzing the wavelet entropy time evolution. Likewise, the dynamic speckle 
phenomenon mentioned above has also been analyzed by using these wavelet 
based entropy concepts [22]. In a recent paper it is investigated the relation 
existing between these two quantiflers — H and NTWS — when they are used 
for analyzing fBm [23]. Figure 1 (top and bottom) shows the mBm and his 
corresponding NTWS, where the same temporal windows, mother wavelet and 
resolution levels were used. 
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2 EXPERIMENTAL SETUP AND DATA ADQUISITION 



The experimental measures were performed in a laboratory by producing ther- 
mal convective turbulence with two electrical heaters in a row. Three different 
turbulence intensities were generated changing the amount of heat dissipated 
for each electrical heater: normal, soft and hard turbulence. Along the laser 
path three electronic thermometers sense the air temperature — T1,T2 and 
T3, see Fig. 2 (bottom). A time series corresponding to the fluctuations of 
the centroid position of a laser beam's spot (wandering) over a screen, af- 
ter propagation through this time-changing laboratory-generated turbulence, 
were recorded with a position sensitive detector located as screen at the end 
of the path. This record consists of 2, 500, 000 spot beam centroid coordinates 
measurements with 500, 000 data for each laboratory-generated turbulence 
condition. Further details of the experiment can be found at Ref. [1]. The 
temperature and signal records can be observed in Figure 2. There, it can be 
observed that the turbulence is increased, and subsequently decreased up to 
recover the initial situation. 



3 RESULTS AND CONCLUSIONS 



In the present work, we employ orthogonal cubic spline functions as mother 
wavelets. Among several alternatives, cubic spline functions are symmetric and 
combine in a suitable proportion smoothness with numerical advantages. They 
have become a recommendable tool for representing natural signals [24, 25]. 
The signal was divided into 606 non-overlapping temporal windows of 4096 
data points. Resolution levels between j = —7 and j = —3 were used to 
calculate both quantifiers. The first two levels (j = —1 and j = —2) were 
dropped to reduce the noise introduced by the sistem, while the lower levels 
were excluded to reduce nonstationary effects as commented in Ref. [1]. 

Figure 3 shows the quantifiers temporal evolution. Both quantifiers reveal that 
when the turbulence is normal the detector can not be able to resolve position 
differences, and electronic noise associated to the detector is observed. The 
NTWS is near one as is expected for a signal generated by a totally random 
process and the a value matches with a white noise. When the turbulence is 
increased the system changes in an abrupt way — see coordinates' graphs at 
Fig. 2. It is interesting to observe the transition between the different intensi- 
ties of turbulence for the signal and its corresponding quantifiers. The Hurst 
parameter discriminates between the other two increased turbulences. It is 
possible, in average, to associate a value ax = 1.17, ay = 1.04 for the soft case 
and ax = 1.62, ay = 1.51 for the hard turbulence. It should be noted that the 
signal has more regularity for the strongest turbulence. The NTWS diminishes 
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notably showing an increment in the order of the system but it is unable to 
distinguish between soft and hard turbulence giving values of NTWS^; = 0.63, 
NTWSy = 0.66 for the soft turbulence and NTWS^ = 0.57, NTWSj, = 0.63 
for the hard turbulence. 

It can be followed by comparing Fig. 2 and Fig. 3 that the behavior of the 
signal is different for both coordinates. Nevertheless, the temporal evolution of 
the quantifiers is very similar. It is also observed that the system has hysteresis 
effect (see Fig. 2 and Fig. 3) as it was expected. 

The mBm model is justified for modeling the dynamics associated to these 
processes. We conclude that the associated scaling exponent changes contin- 
uously with the turbulence strength. In the future a new generalization will 
be considered and studied, the generalized multifractional Brownian motion 
(gmBm) [26] . That processes consider that scaling exponent variations may be 
very erratic and not necessarily a continuous function of time. The latter con- 
dition is a strong limitation in turbulence studies where the scaling exponent 
can change widely from time to time. 
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Fig. 3. Quantifiers temporal evolution: a (top) and NTWS (bottom) for the x (left) 
and y (right) coordinates. 



